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Abstract 

We introduce a family of probabilistic scale-invariant Leibniz-like pyramids and {d + 1)- 
dimensional hyperpyramids (d = 1,2,3,...), with d = 1 corresponding to triangles, d = 2 to 
(tetrahedral) pyramids, and so on. For all values of d, they are characterized by a parameter 
z/ > 0, whose value determines the degree of correlation between N (d + l)-valued random vari- 
ables (d = 1 corresponds to binary variables, ti = 2 to ternary variables, and so on). There are 
{d + 1)^ different events, and the limit z/ — t- oo corresponds to independent random variables, 
in which case each event has a probability l/{d + 1) to occur. The sums of these N {d + 1)- 
valued random variables correspond to a d— dimensional probabilistic model, and generalizes a 
recently proposed one-dimensional (d = 1) model having g— Gaussians (with q = {ly — 2)/(z^ — 1) 
for ly € [1, oo)) as A^ — 7- cx) limit probability distributions for the sum of the N binary variables 
[A. Rodriguez et al, J. Stat. Mech. (2008) P09006; R. Hanel et al, Eur. Phys. J. B 72, 263 
(2009)]. In the z/ — )■ oo limit the d— dimensional multinomial distribution is recovered for the sums, 
which approach a d— dimensional Gaussian distribution for N — ^ oo. For any ly, the conditional 
distributions of the d— dimensional model are shown to yield the corresponding joint distribution 
of the (d — l)-dimensional model with the same i'. For the d = 2 case, we study the joint probabil- 
ity distribution, and identify two classes of marginal distributions, one of them being asymmetric 
and scale-invariant, while the other one is symmetric and only asymptotically scale- invariant. The 
present probabilistic model is proposed as a testing ground for a deeper understanding of the 
necessary and sufficient conditions for having g-Gaussian attractors in the A^ — )■ oo limit, the ul- 
timate goal being a neat mathematical view of the causes clarifying the ubiquitous emergence of 
(/-statistics verified in many natural, artificial and social systems. 

PACS numbers: 05.20.-y,02.50.Cw,05.70.-a 



I. INTRODUCTION 

In a probabilistic context, scale invariance is said to occur when for a set of A^ random 
variables, ^i, ^2, • • • , ^n, with joint probability distribution Pn{^i, ^2, ■ ■ ■ , ^n), the functional 
form of the marginal probabilities of a (A^ — l)-variables subset coincides with its joint 
{N — l)-variables probability distribution, i.e, when 

/ PN{C.l,C.2,---,CN-l,C.N)dC,N =PN-i{Ci,C2,---,Cn-i)- (l) 

In the absence of independence, this relation — which is trivially fulfilled in the case of in- 



'_(jlobal correlations, which is precisely 
5| comes to play an important role. 



dependent random variables — involves the presence o: 
the scenario where nonextensive statistical mechanics 

This theory, also referred to as g— statistics, generalizes the standard Boltzmann-Gibbs 
statistical mechanics, which in turn is appropriate to describe systems which typically present 
local correlations, if any. In this case, the standard Central Limit Theorem (CLT) ensures 
the appearance of Gaussians as attractors in the thermodynamic limit for the sums of 
independent or weakly correlated random variables with finite variance. 

Within the framework of g— statistics, an extension of the CLT, the so called g— generalized 
Central Limit Theorem (g— CLT), has been recently proved (6|,[7| for the case of q— independence 
— a specific class of global correlations — which states that in this case the attractor distri- 
butions in the thermodynamic limit are the so called q— Gaussians, which in d dimensions 



have the form 
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G',(x)=a.e-^-^"^^ xen, (2) 



where x = (xi, . . . , Xd), g is a real parameter with g<l + |,/3isa positive constant, S is a 
positive definite matrix, e^ = [1 + (1 — q)zY^^^^'^'' {e\ = e^), is the so called, q— exponential 
function, whose support is Vtq = M"^ for g > 1 while Qg = {x / x^Y^x < a,^_ A for g < 1, 
and C'^\ = f^ e~^^ ^^dxi ■ ■ ■ dxd, is the normalization constant. The m—th order moments 
of distribution ([2]) are defined for all m if g < 1, and only if g < 1 + ^^^ for g > 1, with 
E[X] = and a covariance matrix given by 

The Gaussians (as well as the independence and the standard CLT) are recovered from 
Eq. (|2]) for g = 1. Another particular instance of g— Gaussian is the uniform distribution, 
which emerges in the g — ?► — 00 limit. 



A number of recent works address the possible relationship between sea' 



relations and g— Gaussians attractors 10Nl4|. Though some discrete 



e invariant cor- 
1 or continuous 



ll| scale-invariant systems have been shown not to have o— Gaussians, but remarkably 
close functions instead, as limiting probability distributions [l2|, the one-dimensional scale- 
invariant model introduced in [131] has been analiticaly shown to yield g— Gaussians in the 
thermodynamic limit. Our goal in the present paper is to study a natural generalization of 
the aforementioned model to higher dimensions. 

In Sec. ini we go through the detailed description of scale invariance. In Sec. [1111 we 
briefly review the one- dimensional model based on Leibniz-like triangles. In Sec. HV] we 
introduce Leibniz-like pyramids and deal with the two-dimensional model. We then explore 
the conditional, marginal and joint probability distributions of the two-dimensional model in 
Sec. lYl We generalize the model to arbitrary dimension in Sec. I VII Finally, we summarize 
our conclusions in Sec. IVIII Some lengthy calculations are developed in the Appendix. 

II. SCALE INVARIANCE 

In order to illustrate the concept of scale invariance, let us consider a statistical model 
consisting of a set of A^ identical and exchangeable binary random variables ^i,^, . . . ,C,n, 
taking values ^^^' and ^^"^^ For the trivial A^ = 1 case there are only 2^ = 2 events in the 
sample space and we have the probability distribution 

Pi (6) = ri,o5(ei - e^) + r,,Mi - e^'^) (4) 

where ri^„ (with n = 0,1, and ri^o + '"o,! = 1) stands for the probability that the variable ^i 
take n times the value $,^'^\ 

For the N = 2 case we have 2^ = 4 different events and the corresponding probability 
distribution 

P2(6, 6) = ^2,0^(^1 - e^' V(6 - e^^^) + r2,Mi - e^'Vfe - e^'^) 

where now the coefficients r2^„ (with n = 0, 1,2, and r2fi + 2r2,i -|- ^2,2 = 1) stand for the 
probability that the value ^^^^ appears n times in the pair (^1,^2)- Note that in order to 
have exchangeable variables the probabilities associated to events {^^'^\^^^^) and (^^^\^*^^-') 
are necessarily the same. 



Let us integrate now expression ([5]) with respect to variable ^2 in order to obtain the 
corresponding marginal distribution 

Pi(ei) = jp2{^i, 6)d6 = (r2,o + r2,i)S{^i - e^) + (r2,i + r2,2)S{^i - ^^'^) (6) 

In order to Eq. ([T]) be fulfilled, i.e., Pi{C,i) = Pi{^i), it is necessary and sufficient that 
r2,o + ^2,1 = ri^o and r2,i + r2,2 = n^i- 

For the A^ = 3 case, with 2^ = 8 events in the sample space, the probability distribution 
reads 

+ r3,25(6 - e^'V(6 - ^^'^ms - e^')) + r3,2<5(ei - e^'Vte - e^^Vfe - e^'^) 
+ r3,25(6 - e^'Vte - e^ms - e^'^) + r3,35(6 - e^'Vte - e^'Vfe - e^) 

(7) 

where r^^n (with n = 0, 1, 2, 3, and r3^o + 3r3_i + 3r3^2 + '"3,3 = 1) stands for the probability that 
value ^^'^^ appears n times in the event (^1,^,^3). Due to the exchangeability of variables, 
there are only A^+1 = 4 different probabilities out of the 2^ = 8 different events. Integrating 
now with respect to ^3 in (^^ we get the marginal distribution 

^2(6,6) = jp3{^u^2,^3)d^3 = (r3,o + r3,i)5(6 - ^^'^)m - ^^'^) 

+ {rs,, + r;,2m, - e^)6{^2 - ^^'^) 

+ (r3,2 + r3,3)<5(ei-e^'V(6-e^'^) (8) 

Now, for marginal distribution ([8]) to coincide with joint distribution ([5]) the relations be- 
tween probabilities r^^o + r^^i = r2,o, "^3,1 + ^3^2 = ''"2,1 and r3^2 + ^3,3 = ''"2,2 rnust hold. 
In the general case, we have 

N 
P^(6,e2,...,e7v) = 5^r;V,n Yl ^i^l-^P)^i^2-d'^)---S{U-^^N^) (9) 

where superindex C runs over the (^) elements set C^ of n— combinations of the A^ elements 
set Af = {1,2,..., A^} and Q (for i G Af) equals ^*^^^ if subindex i is selected by combination 
C and C,^^^ otherwise. The total number of summands in Eq. ([9]) is thus 2^. 



Integrating out with respect to ^n in & one gets the marginal probabihty distribution 



7V-1 



(10) 
which coincides with the joint probabihty distribution pAr_i(^i, ^2, • • • , ^n-i) in case that 

rN,n + rN,n+l = '"Af-l.n- (H) 

This relation is known as the Leibniz rule, and reflects scale-invariance of exchangeable 
binary random variables. We will show a family of models satisfying relation (ITT]) in Sec. IIIII 

Let us consider now a set of A^ identical exchangeable ternary variables, taking values 
^(1)^ ^(2)^ g^]^^ ^(3)_ Pqj, ^Yie N = 1 case we have 3^ = 3 different events and the probability 
distribution 

Pi (6) = n,o,o<5(ei - e^'^) + ri,i,o5(ei - e^'^) + ri,o,Mi - ^^'^) (12) 

where ri^n,m, (with n, m = 0, 1; ^ n + m ^ 1, and ri^o,o + ^^1,1,0 + ^1,0,1 = 1) stands for the 
probability of obtaining n times the value ^'•^^ and m times the value ^^^\ 

For the N = 2 case there are 3^ = 9 different events and the joint probability distribution 
reads 

P2(ei, 6) = ^2,0,0^(^1 - e'^)m - e'^) + r2,i,o5(ei - e^'Vfe - ^^'^) + r2,i,o5(6 - ^^'^)m - e^] 

+ r2,o,Mi - e^'^)^(6 - e^'^) + ^2,0,1^(6 - e^'^)5(6 - e^) + r^xMi - e^'^)5(6 - e^')) 

+ r2,i,i<5(6 - ^^'^m2 - e'^) + r2,2,o'5(6 - e'^m2 - e^'^) + r2,oMi - e^'^)5(6 - e^'^) 

(13) 

where r2,„,m (with n, m = 0, 1, 2; < ri + m ^ 2, and r2,o,o + ^^2,2,0 + ^^2,0,2 + 2(r2,o,i + ^2,1,0 + 
^2,1,1) = 1) stands for the probability of obtaining n times the value ^^'^' and m times the 
value ^^^^ in event (^1,^2)- Due to exchangeability, the number of different probabilities 



reduces to '-^ = 6. 



(Af+l)(Af+2) 
2 

Integrating now with respect to ^2 one gets the marginal distribution 

Pi(6) = / ^'2(6, 6)d6 = (^"2,0,0 + ^^2,0,1 + ^2,1,0)^(6 - ^^^^) 

+ (r2,i,o + r2,i,i + r2,i,i)^(ei-e^'^) 

+ (r2,o,i + r2,o,2 + r2,i,i)6i^i - ^^^^) (14) 



which coincides with the A^ = 1 distribution ( 1T2|) when r2,o,o + ^2,0,1 + ^2,1,0 = ''^1,0,0, ''^2,1,0 + 
^2,1,1 + ^"2,1,1 = ri,i,o and r2,o,i + r2,o,2 + ^2,1,1 = ri,o,i. 
In the general case, we have 

N 
Pn{^1,^2,...,^n)= Y. ^^'".- E ^i^l-^P)^i^2-^P)---SiU-^N^) (15) 

where superindex C runs over the (^^) elements set C^^ of (n,?7i)— combinations of the N 
elements set A/" = {1, 2, ... , A^}, and ^ (for i G Af) equals ^^^•' if subindex i belongs to the 
n elements subset of combination C, ^*^^^ if subindex i belongs to the m elements subset of 
combination C, and $,^^^ otherwise. After integrating with respect to ^at in (TT5|) the following 
marginal probability distribution is obtained 

N-l 

PN~i{(,1, (.2, ■ ■ ■ , ^N-l) = 2_^ {^N,n,m + fN,n+l,m-l + '^AT.m.m-l) 

n+m=0 

X E '^(ei-d'^)5(e2-ef)---'^(e^-i-d'li)- (16) 

'-' ^ ^-^ n , rri 

This distribution coincides with the A^ — 1 variables joint probability distribution if 

^N,n,m + fN,n+l,m-l + fN,n,m-l = fN~l,n,m-l (1 ' J 

We shall refer to condition f lT7|) as the generalized Leibniz rule. We shall show a family of 
models satisfaying such relation in Sec. IIVI 

In the case of A^ identical interchangeable (rf+l) — ary variables taking values ^'^^•*,^*^^-', . . . ,^*-'^^^'', 
the joint probability distribution reads 

AT 
Piv(6, 6, • • • , e^v) = Y. ^N,n,,n,,...,n, J] ^(^1 " ^^(^2 " ^f^) " " " ^i^N ' ^f ) 

(18) 
where superindex C runs over the (^ ^ ) elements set C*^ ^^.....n^ ol ('^i, ''^2, ■ ■ ■ , n^)— combinations 
of the A^ elements set M = {1,2,..., A^} and Q (for i G M) equals ^(•'+^) if subindex i 
belongs to the Uj elements subset corresponding to C, for j = 1, . . . , c?, and ^'^^^ otherwise. 
For c? ^ 2, the joint probabihty distribution (fTSjl satisfy scale invariance condition ([1]) if 
{n = (ni,ra2,...,nrf)): 

^M i^H J L^H -^H fiq^ 

' N,n ^ ' N,n+ei ' ' ' N,n+ed ~ ' N-l,n+ed ^-^^/ 



where e i = e i - 62, Ei = Cj+i - 62, for i = 2, . . . , rf - 1, and Ed = -62; Cj = (0, . . . , 1, . . . , 0) 
being a rf— dimensional vector whose only nonzero component is the i—th one, taking value 
1. Relation fll9p reflects the scale- invariance of exchangeable (rf + l)— nary random variables. 
We shall dedicate the following sections to the detailed description of scale invariant 
probabilistic models. 

III. ONE-DIMENSIONAL MODEL REVISITED 

Let us consider a random experiment which consists in flipping A^ biased coins. We shall 
call X(i) the random variable that counts the number of, say, heads. X(i) G {0, 1, ... , N} 
can be written as 

X(i) =Xi+X2 + --- + Xm (20) 

where each Xj, i = 1, . . . ,N, is a binary random variable taking values 1 (head) with 
probability p, and (tail) with probability (1 — p). As it is well known, considering that 
different throwings are independent of each other, we have Cxx — P{^ ~ P)^iji ^^d the 
probability of obtaining n heads in N trials is given by 

P(X(i) =n)^ p^,„ = (^V"(l - P)^-" (21) 

with X]n=oP^'" ~ -'-' which is no other than the binomial distribution, i.e., X(i) ~ B{N,p), 
with (X(i)) = Np and aj^ = Np{l — p), where the binomial coefficients (^) stand for 
the different ways in which the n heads can be obtained. In other words, though there 
are Q{N) = 2^ different events in the sample space, only X + 1 among them are assigned 
different probability values, namely rj^^n = p"(l — p)^~", n = 0,1, ... ,N, (it is irrelevant 
which specific coins yield head, X(i) only counts the number of them). These selected 
probability values may be displayed in a triangle in the form 

fN,n 



(X = 0) 1 

(X = l) 1-J9 p 

(X = 2) {I-pY {'^-P)P P^ 

(X = 3) {l-pf {^-Pfp i^-p)p'^ P^ 



where the A^-th row displays the A^ + 1 probabihties, r^-n-, for n = 0, 1, . . . , A^. In order to 
get the actual probabilities fl2T]) the above triangle has to be multiplied, element by element, 
by the Pascal triangle. 

Let us now show explicitly the scale invariant character (which otherwise follows trivially 
due to the independent character of the variables) of our probabilistic model. Any of the 
binary variables follows a Bernoulli distribution, pi{Xi = Xj) = p^'(l —pY~^^ with Xj = or 
1, and i = 1, . . . ,N . Thus, as the A^ variables are independent, we have 

N N 

N X^ ^i N -y^^Xj 

P7v(xi,X2,...,X7v) = JJpi(Xi) =p *=1 il-p) '=^ (22) 

Considering now the marginal distribution corresponding to the first (A^ — 1) variables one 

gets 

N-l N-l N-1 / N-1 

1 J2^i N-j2^^ i + J2^' N-ii + j2^i 



xiv=0 



Af-1 AT-l 



p ^=1 (l-p) ^=1 (i + Y^j (23) 

Af-1 7V-1 

^ Xi N -1-^Xi 

p i=l (1-p) i=i =Pn_i{xi,X2,...,XN-i), 



SO a discrete version of Eq. ([T]) is fulfilled and the model is scale-invariant. In addition, it 
is readily seen that coefficients rjsr^n displayed in the above triangle follow the Leibniz rule 
(Eq. (II ip . which, as stated in Sec. HIl serves as an alternative description of scale invariance 
for binary variables). Thus, following Eq. (TTTj) . the sum of two consecutive coefficients in 
any row of the triangle yields the coefficient on top of them. 

Let us recall that the CLT states that, after properly centering and rescaling, one obtains 
a Gaussian (which can be seen as a g— Gaussian with q = 1) distribution out of (l^T]) . namely 
— — ^' A/'(0, 1) for A^ — )■ oo. We shall try now to modify the model so as to obtain 



g— Gaussians with g 7^ 1 in the thermodjTiamic limit 



A. Scale- invariant triangles 

So far, we have considered independent variables. We shall now introduce correlations 
in the model in such a way that the Leibniz triangle rule f lTT]) is preserved — that is, scale 
invariant correlations — by substituting the probabilities rj^^n = p"(l — p)^~"- in i^T\i by 
appropriate ones. In other words, we shall change the above triangle by another one also 
satisfying f lTT]) . As a first attempt, we may resort to the so called Leibniz harmonic triangle 



15|, whose coefficients, defined us 



^N,n^J];^^YuNT = BiN-n + l,n + l); n = 0,l,...,N (24) 

where B{x,y) is the Beta function, may be displayed in triangular form us 



{N = 0) 1 

1 1 



(iV=l) 
(iV = 2) 
(iV = 3) 
(iV = 4) 
(iV = 5) 



1 ^^ ^ ^^ ^ 

^ 4 ^ 12 ^ 12 ^ 4 ^ 

^ 5 ^ 20 ^ 30 ^ 20 ^ 5 ^ 

6 30 60 60 30 6 



Making use of the properties of the Beta function it is a simple task checking that the Leibniz 
triangle above satisfies the rule (TTT!) which, for this reason, is named Leibniz rule. 

Substituting now the rN,n coefficients by the Leibniz coefficients (^^ in Eq. fl2T]) we 
obtain 

with ^n=oPNn = 1' ^hat is, a uniform distribution, which, in the continuum limit, corre- 
sponds, as already mentioned in Sec. HI to a g— Gaussian with q — )■ — oo. 

A family of scale-invariant triangles can be now obtained as properly normalized subtri- 
angles of the Leibniz triangle in the following fashion. Take the central coefficient of any 
even row of the Leibniz triangle and divide the whole triangle by it so as to turn the said co- 
efficient to one. Now take this coefficient as the vertex of a new triangle starting downwards 

10 



from it. The coefficients rjy^ of this new triangle starting from the 2(z/ — 1)— th row of the 
Leibniz triangle can be then expressed as a function of the coefficients rjy^ of the Leibniz 
triangle as 

M _ ' N+2{u-l),n+u-l , . 

' N,n — (T) V^OJ 

'"2(!/-l),!/-l 

As the only transformation we have made on the coefficients of the Leibniz triangle is a 
multiplication by a constant factor, Eq. ( iTTl) still holds for the coefficients (l26l) . thus 



r 



M , ^M _^M 



N,n + ^7V,n+l ~ ''^N-l,n (2' ) 

for any positive integer u. In virtue of Eq. ( l27j) . different coefficients of the triangle are not 

independent of each other and it suffices to specify one element of each row (for instance 

the left side of the triangle) to completely determine the whole triangle. 

Expressing now the Leibniz triangle coefficients in fl26|) in terms of the Beta function as 

in (l24l) . one easily gets 

(,) _ B{N-n + i.,n + i.) _ 

'^^'"^ B{u,u) ' "^^^ ^^^^ 

where now u may take any positive value (see details in |14|). 

Apart from its symmetry, it is worth noticing another remarkable property of coefficientes 
( 128|) . namely, scale invariance condition (127|) guarantees that the set of corresponding prob- 
ability distributions 

&=('')rZ (29) 



n 



associated to the new set of variables 



Xl^] = Xi"^ + Xf ) + ■ ■ ■ + X^^^ (30) 

are well defined, that is, the normalization condition, ^n=oPNn — 1' holds for any value 
of u. It may be shown that (X^^^) = y for all u and a^ (^j = ^^t^^- Concerning the 
binary variables in the sum ( l30l) (taking values and 1 as in the independent case), it may 
be shown that now {X^'"'^) = rj^d = | Vz/, aL,, = | Vz/ and ^L.) (.) = 4"^!, - i = j^r^ 

i i j 

for i ^ j. The constant value of the correlations, which are independent of the system size 
N and the pair of chosen variables X^ , X , reveals the scale- invariant character of the 
model. 

Probability distributions (1291) are shown in Fig. [1] for X = 100 and u = |, 1, |, 2 and |. 

11 



B. Boltzmann and thermodynamic limits 

The limits z/ — )► oo (that we shall call Boltzmann limit for a reason to be clear soon) and 
N ^ oo (thermodynamic limit) are interchangeable in (12 8 p (equivalently in (12^ ). In effect, 
in the first case on has limj^^oo r^n = -^^ ^o 



^^A=(Z]^ (31) 



uiT^'"'''^ Vn/2^ 

and a binomial distribution with p = \ (fair coin), associated to variable X(i) = limi,_5.oo X(^\ 
is recovered, which is consistent with the fact that lim,^_>oo cr^ ,^) ,^) = 0, thus having inde- 

* J 

pendent coins. Then, as the CLT states, by properly centering and rescaling variable ^(i), 
one gets a Gaussian out of (IHT]) in the thermodynamic limit. 

The limit A^ — ?■ oo is much more subtle. It was studied in 13| and later extended in |l4| . 



Depending on the change of variable used when passing from the discrete to the continuous 
model, two families of g— Gaussians are obtained out of fl29l) in the thermodynamic limit, 
with values of q given by q^y = ^j^ for v > 1 oi q^ = |^^ for i/ > 0. Note that two different 

In any case, one gets lim,. 



'-u^oo Hi' 



conjugated g— Gaussians exist for each u > 1 [14 1. 
lim^^ocQu = 1, so the ordinary Gaussian is recovered again. 

To further support our claim, we typify variable X(i) by making the change 

n-^u = ^ (32) 

cr,. 



with a^ = J 4(x+2^) s-iid compare its probability distribution with the corresponding one- 
dimensional gi,— Gaussian ([2]) with q^ = ^5f ^^^ P = g_3 so as to obtain unit variance (see 
Eq. ([3])). Fig. [2] shows (J^p}^[^ versus (n — N/2)/(7,y for z/ = 3 and N = 100 (dots) compared 
with the corresponding g— Gaussian Gq{x) with g = | and unit variance (solid line). The 
agreement is surprisingly good, even for smaller (not shown) values of A^. 



IV. TWO-DIMENSIONAL MODEL 

Let us now generalize our random experiment to one with three different outcomes, which 
is equivalent to throw a biased three-sided dice. Let us label the sides A, B and C and 
suppose that the associated probabilities are P{A) = pi, P{B) = p2, with pi + P2 < 1, and 

12 



P{C) = 1 — Pi — P2- To properly describe our model we need to define a two components 
random variable 

X(2) =X^+X2 + ■■■ + XN (33) 

with X(2) = {X,Y), as a sum of A^ ternary variables Xi = {Xi,Yi), for i = 1,...,N, 
with values Xi = (1,0) (side A), Xi = (0,1) (side B) and Xi = (0,0) (side C). Thus 
X G {0, 1, . . . , A^} counts the number of A's and Y E {0,1,..., N}, with ^ X + Y ^ N, 
counts the number of 5's out of A^ throwings. For independent dices the Xi variables are 
independent of each other, though variables Xi, Yi within each pair are not {{Xi) = {Xf) = 
Pi, (Yi) = (r.2) =p2, o\^x^ = p^{\ - p^)bij, o\^y. =P2{l-P2)Sij, (Tx^Y, = -PiP2Sij), neither 
the variables X and Y, for which (X) = N{pi,p2), and the covariance matrix is given by 

E,., = ivf'"<'"^'' '""" ] (34) 

y -P1P2 P2{l-P2) J 

The probability of having n A's and m B^s, with 0^n + m^N,mN throwings is given 
by 

P{X = n,Y = m)= pM,n,m ={ ^ V>"(1 - Pi - P2)^-"-™ (35) 

\n,mj 

with Y^Q^n+m^NPN,n,m = 1, which is the trinomial distribution, i.e., X(2) ~ T{N,pi,p2), 

where the trinomial coefficients, { ) = f ) f ~") , take into account the degeneracy since 

' \n,m/ \n/ \ m / ^ ° -^ 

among the Vl{N) = 3^ different events in the sample space only '-^ of them have 

different probabilities, namely r]\f^n,m = P"P™(1 ~~ Pi — ^2)^ """^- These probability values 
may be displayed in a pyramid 

(N = 0) 1 



(N = l) 



I-P1-P2 

V2 VI 



(l-pi -P2)^ 
(A^ = 2) P2(1-P1-P2) Pl(l-Pl-P2) 

pI P2P1 pI 

P2{l-Pl-P2)^ Pl{l-Pl-P2)^ 

(AT = 3) P2(1-P1-P2) P2Pl(l-pi -P2) P?(1-P1-P2) 

pI P2P1 P2P1 Pi 
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where each layer, corresponding to a different value of A^, is a triangle with ^ ~^ '^ "*" ^ 
elements, whose rows, from top to bottom, correspond to increasing values oi n + m from 
to A^, while the ascending (descending) diagonals, from left to right (right to left), correspond 
to increasing values of n (m) from to A^. 

In order to get the actual probabilities ( l35l) the above pyramid has to be multiplied by 
the so called Pascal pyramid 



( ) 



iN = 0) 



1 

4 4 



{N = 1) (A^ = 4) 6 12 6 



1 1 

1 

(AT = 2) 2 2 

12 1 



1 



4 12 12 4 
14 6 4 1 

1 

5 5 

10 20 10 

{N = 5) 10 30 30 10 



3 3 

(A^ = 3) 5 20 30 20 5 



3 6 3 
13 3 1 



15 10 10 5 1 



which displays the trinomial coefficients (^ ^) (for the values of n and m within each layer 
the same comments as above hold). 

The joint probability distribution for the independent ternary variables Xj is now 
Pn{xi,X2,...,Xn) = llf=iPi{xi), with pi(xi) = p'^'gy^il - pi - pa)^"""'"^', and Xi = 
(1,0), (0,1) or (0,0). In an analogous fashion as in Eq. (1231) . it may be shown that 
Sx PN(yXi,X2, ■ ■ ■ ,xn) = Pn~i{xi,X2, ■ . ■ ,xn-i), SO the model, as mandated by the in- 
dependence of the variables, still remains scale-invariant. In turn, it is easily seen that 
coefficientes rN^n,m follow the generalized Leibniz rule ( !T7|) . which, as stated in Sec. [TTl 
governs scale invariance for ternary variables. Thus, following Eq. (IT7|) . adding up three 
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elements forming specific triangles within a layer (only those triangles formed by two con- 
secutive elements of any row and the one on top of them), one gets the coefficient just on 
top of them (belonging to the above layer). Fig. [3] shows the first four layers of the pyramid. 
Only the shaded triangles contribute to the generalized Leibniz rule. Two particular sets of 
four elements satisfying the generalized Leibniz rule are shown. 

In analogy with the one-dimensional case, the corresponding bidimensional version of the 
CLT states that the trinomial distribution yields a bidimensional Gaussian in the thermo- 
dynamic limit, that is, for A^ — )■ cxo we have ^^' _ ^^^ ~ A/'(0, ;^S), with S already given 
in f p^ . Our next step is to analyze under what conditions this type of model can provide 
bidimensional g— Gaussians. 



A. Scale- invariant pyramids 



As we did in Sec. IIIIAl we shall now introduce scale-invariant correlations in the model, 
so different dice throwings will no longer be independent. With this purpose, we may extend 
the Leibniz triangle to a Leibniz-like (tetrahedral) pyramid: 



N,n,m 



(iV = 0) 



(iV = l) 



(iV = 2) 



iN = 3) 



1 

3 

1 1 
3 3 



12 



J_ 
12 



J_ 
12 



10 

30 30 
1 1 



30 



60 



30 



J_ 
10 



30 



J_ 
30 



10 



(AT = 4) 



(AT = 5) 











1 
15 
















1 
60 




1 
60 












1 
90 




1 
180 




1 
90 








1 

60 




1 

180 




1 

180 




1 
60 






1 


1 




1 




1 


1 






15 


60 

1 
210 


1 
105 


90 

1 

21 

1 
420 


1 
105 


60 

1 
210 


15 






1 
210 




1 
630 




1 
630 




1 
210 
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1 




1 
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1 






105 


420 




630 




420 


105 
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1 




1 




1 




1 


1 



21 



105 



210 



210 



105 



21 
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The coefficients are given by 

^•"'- (Ar + 2)(iV + l)(„g- ^ ' 

They satisfy the generahzed Leibniz condition flT7|) . and are defined in such a way that, when 
multiplying element by element by the Pascal pyramid, we get for the actual probabilities a 
uniform distribution: 

PN,n,m- y^^^J N,n,m (^N + 2){N+1) ^^ 

with J2n+m=oPNnm — 1; which again leads as to a (bidimensional) g— Gaussian with q — )■ 
— oo. 

As we did with the triangles, we can now get a family of scale-invariant pyramids out 
of the Leibniz-like pyramid. As the only layers of the pyramid with a central coefficient 
are those multiple of 3, we shall descend 3 by 3 layers from the top of the Leibniz pyramid 
and divide the whole pyramid by the corresponding central element, that is, the element 
''^3(1-1) u~iu-i '^^ ^^^ "^(^ ~ l)~th layer. This coefficient so turns to unity and we can start a 
new pyramid downwards from it, with coefficients given by 

(y) A''+3(!^— l),n+j/— l,m+i/— 1 ('i<;i\ 

''"N,n,m — (1) ■ W°J 

By construction, the family of pyramids f l38|) satisfy the generalized Leibniz rule f lT7|l (see 
comment following Eq. ( l26ll ). hence 

' N,n,m ^ ' N,n+l,m-l ^ ' N,n,m-1 ~ ' N-l,n,m-l V'^^J 

for any positive integer u. Due to restriction ( 139|) on the pyramid coefficients, the whole 
pyramid may be determined by only specifying the elements of one face. 

Expressing the trinomial coefficient in (!36l) as a product of binomial coefficients and 
making use of the property of the Beta function expressed in (El]), the coefficients of the 
Leibniz-like pyramid can alternatively be expressed as 

^^,l,m = ^^^^''^"^^^ ^(A^ -n-m + l,m + l)B{N -n + l,n + l) (40) 

Introducing now (l40l) in (!38l) . after some algebra, one finally gets 



(j,) B{N — n — m + v,n + 171 + 2u)B{n + u,m + u) 

B{iy,u)B{u,2iy) 



(jy) ijiiv — lb — lib -T i^, lb -r lib -r ^iyj±jfib -T 1^, lib -r ly J / a-i\ 

^N'n,m = D/,, ,AD/,. n,.\ 5 '^ > (41) 
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where now v can take on any positive value. 

Multiplying the Pascal pjrraniid by the family of pyramids (HTj) . one gets a family of 
probability distributions 

VN,n,ni = L ^ ^1 ^N,n,m (42) 

with ^n+m=o Pnu m ~ ^ ^'^^ ^^Y ^ i associatcd to the new family of statistical variables 

xJj)=xf) + xf) + ... + x(;) (43) 

Fig. m shows probability distributions ( H2|) for A^ = 50 and z/ = | , 1 (uniform distribution 
in the triangle Q ^n + m ^ N), 2 and 5. 



It may be shown that {xXi) = ("f"; "f") ^.nd the covariance matrix is 



,(,) _ N{N + 3z/) / 2 -1 
9(1 + 3//) _i 2 



^2! = ^i)f^£v^ I : ; I (44) 



Concerning the ternary variables X^ = (X^'^ ,Y^ ) in the sum ( 143|) (which take the same 
values as the ternary independent variables Xi in ( l33l) ). they are no longer independent 
of each other, their first and second order moments being given by {X^ ) = {Y^ ) = 

'1,0,0 — 3 ^'■^^ x('') ~ y.C") ~ 9 ^'^' x'^'x*"^ ~ '^v('')y('') ~ '2,0,0 9 — 9(1+3;^) ^"^ ^ / ■>'' 

i i i j * J 

cr^.y. = —1 and cr^ („) („) = Tg 1 q — h = ~^7tttt ^^^ ^ 7^ .^- '^^^ scale-invariant character of 



the correlations that we have introduce becomes, as before, apparent. 



B. Boltzmann limit 

Taking limits in f l4T]) one gets lim,y_j.oo ^jvn m~ 3^? so in the Boltzmann limit the trinomial 
distribution 

withpi = P2 = 1— Pi— P2 = I is obtained (as well as covariance matrix (H4|) . in the Boltzmann 
limit, coincides with covariance matrix flMj) for the same values of p and q), which is again 

consistent with the fact that limi^^oo ^l(>') v^{-) = limi^-^oo (^yMyi-^) = lii^i^^oo '^l(>')y(w = 0, so 
independent dices are recovered. Applying now the CLT, a Gaussian distribution appears 
again in the thermodynamic limit. 

We will devote Sec. |V] to the study of probability distribution ( H2|) . in particular its 
thermodynamic limit, in order to check if it is a bidimensional g— Gaussian. 
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C. Entropy 



As it was shown in Ref. 13|], the family of distributions fl29|) based on scale-invariant 



triangles (1281) . are properly described by the Boltzmann-Gibbs entropy. It has been re- 



cently shown (see 16|] and references therein) that this is so due to the fact the number of 
microstates increases exponentially with the system size (fi(A^) = 2^). In turn, the non- 
additive entropy Sq with g 7^ 1 [1] (the q = 1 case corresponds to the Boltzmann entropy) 
turns out to be appropriate (in the sense that the entropy is extensive) to describe systems 
: or which a relevant fraction of the degrees of freedom vanishes in the thermodynamic limit 

~l n 

Iq . Such is the case for the binary models studied in [17|, corresponding to triangles where 
the nonzero probabilities restrict, in the thermodynamic limit, to a strip of size b on the 
triangle, being well described by the nonadditive entropy with g = 1 — ^ (more precisely, for 
this value of q, Sq is extensive). 

For the ternary model (j^2|) we also have exponential increase of the phase volume {^{N) = 
3^), so, once again, the Boltzmann entropy is expected to be extensive. The nonadditive 
entropy for the family of pyramids f HTl) is given by 

^ - 1 („:j «i.)' 

^(.) ^ n+m=0 ^^g^ 

for g ^ 1 and Sj''^ = -J2n+m=o (i)^/v,l,mln (rJ^J^.^^j for g = 1. Fig. [5] shows g-entropy 
046 p as a function of A^ for different values of z/ and g. As predicted, independently of the 
value of z/, the value of g which makes the entropy Sg extensive is gent = 1- 

V. JOINT, CONDITIONAL AND MARGINAL DISTRIBUTIONS 

In order to get a deeper insight into the two-dimensional probability distribution f l42|) we 
shall start by studying its associated one-dimensional marginal distributions since, as stated 
in Ref. 8(], in case these distributinos were g— Gaussians, the two-dimensional one would 

also be g— Gaussian, with values of g related by j^ — = j^ h 1, where gi^ (g2d) stands for 

the value of g corresponding to the one- (two-) dimensional g— Gaussian. 

By construction, probability distribution (H2|) inherits the highly symmetrical charac- 
ter of the trinomial coefficients (any layer of the Pascal pyramid is a triangle symmetric 
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with respect to any of its heights) and so p%]^^^ = p%]^,^ = p^^^^-n-m^m = P%]n,N-n-m = 
PNm N-n-m — Pnu N-n-m (^^^ Same property holds for coefficients (HT]) ). Thus the marginal 

1- , -1 ,• ~M sr^N—n (u) ~(u) v-^TV— m (u) j ~(u) v-^ (u) 

distributions pV,„ = Em=oPkn,m' PnU = E„=0 PN'n,m and p'^'i = Y.n+m=lPN,n,m^ ^SSOCl- 

ated to variables n, m and I = n + m E {0, 1, . . . , A^} (which, as shown in Fig. El correspond 
respectively to adding up probabilities along the ascending diagonals, descending diagonals 
or rows of the layers of the pyramid) are identical. Fig. [7] shows probability distributions 
p}^[^ for A^ = 100 and z/ = l/2, 1, |, 2, | and 5. The first apparent property of these dis- 
tributions is its nonsymmetric character, their first moments being given by (X) = y and, 
as easily deduced from (jH]), aj^ = gn^J) ■ This asymmetry vanishes when z/ — > oo (as 
stated in Sec. IIVBI in the Boltzmann limit, bidimensional distribution ( l42l) tends to a sym- 
metric trinomial distribution, thus having symmetric binomial marginal distributions) but 
remains when increasing A^ (it may be shown that the skewness coefficient assymptotically 
approaches a finite value in the thermodynamic limit). 

A new family of nonsymmetric triangles can be defined from the marginal distributions 
p]^^ in the form 



~{u) _ PN,n 



a 



(47) 



The V = 1 instance of triangles (H7j) is 



f(i) 

' N,n 







2 
3 




1 
3 












1 
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1 
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2 
5 
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10 
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1 
10 






1 
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4 
15 




1 
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2 
15 




1 
15 




2 5 




4 




1 




2 




1 


7 21 




21 




7 




21 




21 



Though asymmetric, triangles ( HTj) fulfill scale-invariance condition ( l271l . In fact, they can 
alternatively be obtained applying rule ( 127]) to the edge rjy q q = f}^\j of the corresponding 
pyramid. 

In addition, for integer values of z/, triangles Wf\ may be obtained as nonsymmetric, 
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properly rescaled subtriangles of the family (128|) in the form 



41 = "^T (48) 



' u,0 



From ( l48l) . the marginal probabilities pAr_„ can be related to the one-dimensional family ( !29|l 
in the form 

H ^ 1 (iV-n + l)---(iV-n + z.) (.) 



which, for 1^ = 1, reduces to the straight line 



'^"?"=(ivTT)(ivT2)('^-" + ^) ("' 



shown in Fig. [7l 

The nonsymmetric character of marginal distributions p]^^ makes them not good candi- 
ates to yield g— Gaussians in the thermodynamic limit. Nevertheless, there is another yet 
unexplored direction for the marginal distributions. We shall now add up coefficients in 
the vertical direction of each layer, which corresponds to calculate the marginal distribution 
PN,k = Yjn-m=kPN,n,m associatcd to the variable k = n - m e {-N, . . . ,0,. . . ,N} (see 
Fig. [6]). Fig. [8] shows marginal distributions p}^\. for A^ = 100 and z/ = i, 1, |, 2 and |. 
Contrary to the former case, these new marginal distributions are symmetric, due to the 
aforementioned fact that the layers of the pyramid are symmetric triangles with respect to 
their heights. From these distributions and dividing by the appropiate binomial coefficient, 
we can now define a new family of triangles in the form 



f'2N,N+k = / 27V \ ' k = -N,. . . ,0,. . .,N (51) 

\N+k) 



with only even labelled rows. The z/ = 2 instance of triangles (|5T|) is 
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^(2) 

' 2N,N+k 
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Though symmetric triangles (151 p share their sides with one of the sides of triangles f l47p 
and the edge of pyramids (HTI) . i.e., r^^^ = fjy-^ = rjyQ q, it should be stressed that they are 
not scale invariant. To properly define scale invariance in triangles f l5T]) . Leibniz rule f lTTj) 
can recursively be applied to connect coefficients separated two rows, thus obtaining r^^ ^ + 
2^27Vfc+i + ^2Nk+2 — ^2N~2k- This relation is not fulfilled by triangles (ISTI) . Nevertheless, it 
is asymptotically fulfilled in the thermodynamic limit. In effect, we shall define the quotient 

' 2N,k ' ^' 2Ar,fc+l "•" ' 2Af,fc+2 

and show that MmN^oo Pn n = 1- FiS- E] shows quotient ( l52l) versus k/N for i/ = 2 and 
A^ = 50, 100 and 200. An oscillating trend around the value 1 is observed, the quotient 

(2) 

being closer to 1 for the central values of k. Increasing A^ makes pj^\. closer to 1. We then 
verify that triangles flSTl) are asymptotically scale invariant. 

There is, however, an alternative way of obtaining the same above distributions p% f.- 
Let us think of our ternary variables as being scalars — instead of vectors, as imposed 
in (H3l) — having values —1, and 1. We shall denote them by ^^ , i = 1,...,N, and 



consider the sum of them ^^'^' = X]i=iC > which is a scalar variable having values in 
{— A^, . . . , 0, . . . , A^}. If we associate the probabilities of pyramid P2l) to the probability of 
having n variables with value 1 and m variables with value —1, the probability distribution of 
variable ^^''^ exactly coincides with the marginal distribution p}^^^ of distribution ( 142|) . Thus, 
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in our model the dimension associated to the random variables plays a secondary role, the 
specific form of the correlations among them being of major importance. Some numerical 
results seem to indicate that marginal distributions pl,\, do not yield g— Gaussians in the 
thermodynamic limit either. Nevertheless, the difficulty to reach high values of A^ (due 
to the rapidly increasing value of the trinomial coefficients) doesn't allow us to establish a 
definite conclusion. 

Let us turn now to the study of the conditional distributions of probability distribution 
f H2|) . which are defined in the way 

PN,n\m = ^(^;^-^ n = 0,...,N-m (53) 

PN,m 

and an equivalent expression for pj^v«in' which, for symmetry considerations, is identical. 
We show in the Appendix that 

PN,n\m ^ PN-m,n \^^) 

i.e., the conditional distributions of the family of distributions ( 142|) coincide with the family 
(1291) . already shown in Fig. [1], for an appropiate system size. 

Thus, bidimensional distribution ( H2l) contains sections which, as stated in Sec. IIIIBl are 
g— Gaussians in the thermodynamic limit, which is a necessary, but by no means sufficient 
condition for probability distributions fH2|) to be two-dimensional g— Gaussians. 

To further explore the nature of distribution P2|) we will do the following. Let us first 
typify variable X^^^ by making the linear change 

^"\("]^a( "-'"'] (55) 

yv J ym-N/3 J 

with jacobian (det A)~^, where matrix A is defined in such a way that {T,Z!)~^ = A^A, with 




(2)' 
(2) 



S^^] given in (jS]), thus 




Vector random variable U = {u, vY defined in ( l55l) is centered, i.e., E\U] = A{E[X/^}] —jl) = 
0, where jl = (y, y) , and has identity covariance matrix, since E[UU'^] = E[A{X^^!^^ — 
/x)(xg -ilfA^] = AE[(X{2)^ -/i)(xg -/i)^]A^ = AsJ^Ja^ = I. In addition, change §B 
has the extra effect of converting triangle ^ n + m ^ N oi vertices (0, 0), (A^, 0) and (0, A^), 
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into a centered in the origin and equilateral triangle of vertices a(— 1, 0), a(l/2, — -\/3/2) and 
a(l/2, a/3/2), with a = \/ 7v7F+^' i^o^e appropiate for our purpose. 



In the hypothesis that distribution P2|) yields a bidimensional g,y— Gaussian with g^, 



i/-2 

in the thermodynamic limit, we shall compare probability distribution of variable U with 
a typified bidimensional g;^— Gaussian ([2]) with (3 = g_^ and S = J (thus having identity 
covariance matrix as ensured by Eq. ([3])). Fig. ITO] plots {det A)~^p}^'^^ versus U (dots) 
compared with the corresponding bidimensional gi,-Gaussian (solid surface), for A^ = 50 and 
u = 2, q^ = (left), u = 10, q,y = ^ (center) and u = 20, q^ = ^ (right). For z/ = 2 the 
mismatch is evident, the results not improving when increasing A^, neither by changing the 
value of q. Nevertheless, when increasing u the fitting improves, which is due to the fact that, 
as stated in Sec. IIVBI family (H2l) approaches a trinomial distribution in the Boltzmann 
limit, which in turn approaches a Gaussian distribution in the thermodynamic limit, while 
the corresponding g^— Gaussian also approaches a Gaussian since lim,^_>oo lu = ^■ 

Thus, bidimensional g— Gaussians with generic values of g are elusive in the present model. 
A possible reason for it is the inadequacy of linear change (l55l) . Instead, a nonlinear (and 
highly not trivial) change is needed, which transforms the triangle ^ m + n ^ N into a 
circle (either for all values of A^, or at least for increasing A^). Research along this line is in 
progress. 

VI. GENERALIZATION OF THE MODEL TO ARBITRARY DIMENSION 

We may now throw a four-sided dice (i.e., a tetrahedric dice), with associated probabilities 
Pi, P2, Ps (with pi + P2 + P3 < 1), and I — pi — P2 — Ps for the different sides. We need now 
a three components variable 

X(3)=Xi + X2 + --- + X^ (57) 

with X(3) = {X,Y,Z), defined as a sum of A^ quaternary variables X, = {Xi,Yi, Zi), for 
i = 1, . . . ,N, taking the four possible values (1, 0, 0), (0, 1, 0), (0, 0, 1) and (0, 0, 0) associated 
to the different sides, so that X, Y and Z count the number of appearences of three of the 
sides out of A^ throwings. Variable X(3) follows the tetranomial distribution 

P{X = n,Y = m,Z = l)=( ^ X'lp^piil - Pi - P2 - p^)""—^-' (58) 
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where ^ n + m + / ^ iV and the tetranomial coefficients, ( ^ ,) = D f^"") f^"""'") , 

^ ^ ' \n,m,U \n/ \ m / \ I /' 

stand for the number of ways in which the same resuh (X, Y, Z) = {n, m, I) can be obtained 

after A^ throwings. We have now n(N) = 4^ different events, though only '^ ^ '^ - 

different probabihty values, namely rN,n,m,i = PiP^Psll — Pi ^ P2 — ps)^"""'""'. In order 

to display them (what we will not do) we would need to use a g elements 

pyramid for each value of A^ (in the same way as we used a ^ '^ '^ '^ ' elements triangle for 
each value of A^ for the trinomial distribution and a A^ + 1 elements row for each value of 
A^ for the binomial distribution). Putting all the pyramids together makes a hyperpyramid 
(in the same way as we got a pyramid made of triangles for the trinomial distribution and a 
triangle made of rows for the binomial distribution). In order to get the actual probabilities 
(j58|) . this hyperpyramid should be multiplied by what we may call Pascal hyperpyramid 
{N = 0) iN=l) iN = 2) iN = 3) {N = A) 



1 2 

2 2 
1 

11 1 

2 2 
12 1 



4 
1 4 4 

3 6 

3 3 12 12 

6 12 6 
3 

6 6 4 

3 6 3 12 12 

12 24 12 
1 4 12 12 4 

3 3 
3 6 3 1 

13 3 1 4 4 

6 12 6 

4 12 12 4 

14 6 4 1 
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which displays the tetranomial coefficients. Now, layers from top to bottom within A^— th 
pyramid, correspond to constant values oi n + m + I from to A^. In addition, rows from 
bottom to top within each layer correspond to constant values of / from to / + m + n, while 
indexes m and n behave in the same way as in the Leibniz pyramid. 
We may define now the Leibniz hyperpyramid as 



.(1) ^ ^ 

N,n,m,l / j>T i\/Ar oX/Ar i o\ ^ N 



{N + l)iN + 2)iN + 3){JJ 
whose coefficients may be displayed as 

(TV = 0) {N = 1) {N = 2) {N = 3) {N = 4) 



(59) 
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So defined, when multiplying Pascal hyperpyramid by Leibniz hyperpyramid one obtains 
the uniform tridimensional distribution given by 

PN,n,m,i = (Ar + i)(iV + 2)(Ar + 3) ^^°^ 

In adittion, a family of hyperpyramids may be extracted from the Leibniz hyperpyramid 
flM?]) in the way 

{u) _ ' N+4:{u-l),n+u-l,m+u-l,l+u-l .„..s 

'^N,n,mJ — (Ti \^^) 



N,n,m,l (1) 
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^4(l/-l),l/-l,!/-l,l/-l 



for any positive integer u, which satisfy a furtherly generahzed Leibniz rule 



''"N,n,m,l + ''"N,n+l,m-l,l + ''" N,n,m-l,l+l + ''^N,n,m-l,l ~ "^ N-l,n,m-l,l^ 



(62) 



which states that adding up certain four elements subpyramids of the A^— th pyramid (anal- 
ogous to the subpyramids of the Leibniz-like pyramid depicted in Fig. [3]) one gets a corre- 
sponding element of the (A^ — 1)— th pyramid. Relation f l62p is the c? = 3 instance of relation 

In order to get an explicit expression for coefficients f l6T|) . we shall first rewrite Leibniz-like 
triangle coefficients (128|1 in the form 

(,) _ T{2u)V{N-n + v)V{n + v) 



N,n -^21 



l^ 



T{N + 2iy) 

and correspondingly, coefficients ( HTl) of the Leibniz-like pyramids in the form 
(^) _ r(3z/) T{N -n-m + u)T{n + z/)r(m + u) 



(63) 



(64) 



which naturally leads us to rewrite the coefficients fl6Tl) of the family of hyperpyramids as 



M) 

N,n,m,l 'Pi I 



r(4z/) r(A^ -n-m-l + p)V{n + z/)r(m + z/)r(/ + v) 



y 



which satisfy relation fl62|) for any z/ > 0. 

Generalizing former results, the corresponding family of probability distributions 



(65) 



PN,n,mJ 



N 
n, m, I 



S^) 



(66) 



associated to random variable X)^} = ^ .^-^ X^^^^with {X^^l^h = (^, ^, ^) for all u and 



covariance matrix 



'(3) 16(1 + 4z/) 



/ 



v 



3 -1 -l\ 
-1 3 -1 
-1-1 3 / 



(67) 



yields a tetranomial distribution, with pi = p2 = Ps = 1 — Pi — P2 — Ps = |, in the Boltzmann 
limit: 

where independent dices are recovered. 



26 



We are finally prepared to extend the above structure to higher dimensions. To de- 
scribe a random experiment consisting of throwing a (rf + l)-sided dice N times, we need a 
(i— dimensional random variable X(^d) = (-'^i, -^^2, • • • , ^d), with 

X(,)=Xi + X2 + --- + X^ (69) 

— * 

where Xj = {Xij,X2j, . . . , X^j), j = 1, . . . ,N, is the li-dimensional, {d + l)-valued random 
variable associated to the j— th throw, taking values Ci = (1, 0, . . . , 0) for side labelled 1, 
62 = (0, 1, ... , 0) for side labelled 2, . . . , e^ = (0, . . . , 1) for side labelled d, and = (0, . . . , 0) 
for side labelled (d + l) {d = 1 corresponds to binary variables, d = 2 corresponds to ternary 
variables, and so on). Thus, variable Xj = X]i=i -^ij' ^'^^ i = 1, . ■ ■ ,d, counts the number of 
appearances of side labelled i out of N throwings. 
Variable X(^a) follows a multinomial distribution 

P(Xi = n,,X2 = n2,...,Xd = n,)=( ^ ^^^2' " ■■Pd'i^-Pi Pd)'""^--"'* 

\ni,n2,...,ndj 

(70) 
where Pi stands for the probability of obtaining side i, ^ ni + ^2 + ■ ■ ■ + n^ ^ N, 

and the multinomial coefficients, ( ) = — j — j rjir tt, stand for the different 

ways to obtain the result (ni,n2, . . . ,'^d)- Thus, there are only — ^, '"'^ different 

probability values, namely r7v,ni,n2v,nd — pTpT ' ' 'Pd'^i^ — pi — • ■ ■ — pd)^^"^ "'^ among 
the Q{N) = {d + 1)^ events of the sample space. 

The {d + 1) -dimensional Leibniz hyperpyramid is given by 



N,n.,n2,...,n, (iV + 1) (iV + 2) . . . (iV + d) („^^ J... J 



^N^m.no n,i /AT , -i\/Ar , r^\ /AT , 7\ / N \ V'-'-J 



from which the family of hyperpyramids 

(v) _ N+{d+l){v-l),ni+v-l,n2+u-l,...,nci+v-l 

'N,ni,n2,...,nd ~ (i) 

^{d+l){u-l),u-l,i/-l,...,L'-l 



r{{d + l)u) r{N -ni-n2 n^ + u)r{ni + u)r{n2 + u) ■ ■ ■ r(n^ + u) 

T{uy+^ T{N+{d+l)u) 

(72) 



with associated probabilities 



» _ / ^ \J'^) 



PN,nun2,...,n, [^^^^ n2, • • • , uj ^^'"i'"^'-'"'^ ^^^^ 
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for variable X, 



(d) 



Eti^t^ is obtained, with (xgh = (J^ 



covariance matrix 



/ 



'id) 



N{N+{d+l)u) 

(rf+ 1)2(1 + (rf+l)z/) 



'id)' 

d -1 
-1 d 



\d+l- 



-^^ for all z/, and 



A 



\ 



■1 d 



(74) 



/ 



(^) 



(concerning the component variables X^ ■ 



l,...,d. j = l,...,iV, onehas(xf7) 



a 



X, 



(",) ~ (d+l)2 ' 



a 



i.j i,k 



r, 



2,0,. ..,0 



for i 7^ /, and a 



,H 



1 

(rf+TF 
1 



for j 7^ A;, o"^ 



d+i' 
1 



^(^)^(^) — ' 2,1,0,. ..,0 (d+l)2 



i,j fc,i 



d+l)2{l + (d+l)!/)' ^'-'^ J T l^l '^ y('') v('-') (d+l)2 ' 

(dTTFIiWM ' fo^ ^ 7^ ^ and J ^ /.) 



The generalized scale invariance rule, followed by coefficients f[72|) for c? ^ 2, is given in 
Eq. ^. 



VII. CONCLUSIONS 



We have generalized to an arbitrary dimensron a one drnrensronal discrete probabiUstic 
model ffist proposed in [131 which, for one dimension, yields g— Gaussians in the thermody- 
namic limit. These functions, which play a central role in nonextensive statistical mechacics, 
appear in a plethora of applications to natural, social and artificial systems. 

Though we have obtained two-dimensional distributions which contain one-dimensional 
conditional distributions yielding one-dimensional g— Gaussians in the thermodynamic limit, 
our calculations seem to indicate that the model does not yield, for generic g, bidimensional 
g— Gaussians as limiting probability distributions for A^ — )■ oo. In addition, making use of the 
corresponding generalization of Eq. ( 154|) . our model contains one- dimensional g— Gaussians 
as conditional probability distributions for any starting dimension d (corresponding in turn 
to sums of {d + 1)— ary variables). 

We are thus lead to the conclusion that the case of binary variables is special. Indeed, 
it appears to be the only one yielding g— Gaussians in the thermodynamic limit in a simple 
manner. For more complex random variables, an adaptation seems necessary in what con- 
cerns the boundaries that emerge. More precisely, a difficulty exists to match the domain 
where the random variables of our model take their values (e.g., a triangle or a tetrahedron 
for dimensions 2 or 3 respectively) with the support of the g— Gaussians, which is the inte- 
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rior of an ellipse or of an ellipsoid in dimensions 2 or 3 respectively, for the case q < 1; or 



the whole euclidian space for g > 1 (in which case, as shown in Ref. IJ] for dimension 1, 
a complicated change of variables is needed in order to transform the bounded support of 
the discret model into the whole real axis). This difficulty becomes apparent when trying 



to use a sort of generalization of the Laplace-de Finetti theorem [18|, which allows for a 
representation of exchangeable stochastic processes, and has been successfully applied to 
the case of binary random variables [1^ . 

Though the formulation of a simple probabilistic model yielding multidimensional 
g— Gaussians in the thermodynamic limit still remains an open question, we have intro- 
duced a very rich scale-invariant multidimensional model which further support the hy- 
pothesis that scale invariance is possibly a necessary — though definitely not sufficient — 
condition for g— independence, which in turn guarantees the appearence of g— Gaussian 
attractors. 
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APPENDIX 

In order to prove relation fl5^ . we will substitute Eq. fl53|) . with Py^n^ given in fH2|) . in 
the l.h.s. and Eq. (129|) in the r.h.s. of (15^ . to obtain 

^ \ H (N -m\ (^) 



N,n\m N-m / ,, x N-m /,, x \ n T ^-'^." ^ ^ 



^N,n,m 



71=0 ^ ' ^ n=0 ^ ^ 

where we have made use of the relation (^ ^) = (^) ( ""') . Thus, Eq. (15^ is fulfilled 



whenever 



~ ' N-m,n \'^) 



Z^ I ^ r N,n,m. 

n=0 ^ ^ 
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which estabhshes a relationship between properly rescaled elements of the pyramid and the 
elements of the triangle with the same value of v. 

Introducing now in (175]1 expressions (1^ and (IMI) for the coefficients of the triangles and 
pyramids respectively, and eliminating terms which cancel out, one gets 



T{2v) 



N-m\ T{ufT{N + 2u) 



( \T{N -n-m + v 

n=0 \ ^ / 



(77) 



which, making again use of Eq. f l63p . is trivially fulfilled 

r(2z/) ^/AT-m 
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FIG. 1. (Color online) Probability distributions ([29D for N = W0 andu = ^, 1 



|, 2 and |. 
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FIG. 2. cr^p}^^ versus {n — N/2)/ay for v = 3 and N = 100 (dots) compared with the corresponding 
g— Gaussian with q = 2 ^^*^ ^^it variance (sohd hue). 
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3.0.3 




2,0,2 



FIG. 3. Structure of the pyramid formed by coefficients r]\[^n,m = p"'q"^{l — p — qr)^~"~™ (or 
coefficients of the family (j4ip ). Only the shaded triangles contribute to the generalized Leibniz 
rule. The four coefficients involved in Eq. (J17p form in turn a subpyramid of the pyramid. Two of 
such subpyramids are indicated in de Figure. 
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FIG. 4. Probability distributions (^ for A^ = 50 and z^ = ^ (top left), 1 (top right), 2 (bottom 
left) and 5 (bottom right). 



35 




n I ^ I ^ I r 
20 40 60 



1 ' \ ' \ ' r 

20 40 60 80 



FIG. 5. (Color online) Tsallis entropy (gSl) for q =0.8, 0.9, 1, 1.1 and 1.2 of family ([H]) for u = 1, 
2, 3 and 4. For all values of u we have gent = 1- 
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n-m = -2 n-m = n-m = 2 



m = 1 
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n = 1 
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n = 2 



n = 3 
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n + m = 



n-\- m = 1 



n + m = 2 



P3.0.1 



P3.1.0 



K3,0,Z 



C") ,- 

P3,l,l 



C") 
Ps.i.o 



n + m = 3 > P3.0.3 P3.1.2 ' P3.2.1 P3.3.0 



t t 



n — m = —3 n — m = —1 n — m=l n — ni = 3 

FIG. 6. A^ = 3 layer of pyramid ()42p . The three heights of the triangle are shown. As the elements 
of the triangle are symmetric with respect to any of them, there are only three different elements 
in this layer. Left: ascending (descending) diagonals of the layer correspond to constant values of 
n (m). Right: rows (columns) of the layer correspond to constant values oi n + m (n — m). 
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FIG. 7. (Color online) Marginal probability distributions i?)v„ for A^ = 100 and u 
and 5. 
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FIG. 8. (Color online) Marginal probability distributions p}!j \^ for N = 100 and 
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FIG. 9. (Color online) Quotient ([52]) versus k/N ioi v = 2 and N = 50, 100 and 200. When 

(2) 

increasing A^, yO]^^ approaches 1. 
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FIG. 10. (det A) ^p^[^ ^ versus U (dots) compared with the corresponding bidimensional q^- 
Gaussian (solid surface) for A^ = 50 and v = 2, q^ = {) (left), v = 1{), q^ = ^ (center) and v = 20, 
<ly = ^ (right). 
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